Mandelbrot Law of Evolving Networks 
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Degree distributions of many real networks are known to follow the Mandelbrot law, which can be considered 
as an extension of the power law and is determined by not only the power-law exponent, but also the shifting 
coefficient. Although the shifting coefficient highly affects the shape of distribution, it receives less attention in 
the literature and in fact, mainstream analytical method based on backward or forward difference will lead to 
considerable deviations to its value. In this Letter, we show that the degree distribution of a growing network 
with linear preferential attachment approximately follows the Mandelbrot law. We propose an analytical method 
based on a recursive formula that can obtain a more accurate expression of the shifting coefficient. Simulations 
demonstrate the advantages of our method. This work provides a possible mechanism leading to the Mandelbrot 
law of evolving networks, and refines the mainstream analytical methods for the shifting coefficient. 
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Many systems can be described as networks yj-lj], in 
which, the nodes correspond to the elements and the links to 
the relations between elements. Uncovering the mechanisms 
underlying the structural features of real networks is one of 
the most significant challenges in network science. Two pi- 
oneering models, respectively for small-world |5] and scale- 
free networks |6], give explanations for many real phenom- 
ena, such as, the logarithmic growth of average distance, the 
power-law degree distribution, and the high clustering coeffi- 
cient. With the idea of 'rich get richer', the Barabasi- Albert 
(BA) network |6] embodies two mechanisms: growth and 
preferential attachment. That is, at each time step, a new node 
is added and connected to a few old nodes with probability 
proportional to their degree as: 



n(fe,) = fc./^fcj 



(1) 



where ki is the degree of node i, and j runs over all old nodes. 
The analytical solution of the degree distribution. 



p{k) = 2m^k~^, 



(2) 



can be obtained by applying the mean-field approximation [^, 
0], in which 2m is the average degree of the network. 

Unfortunately, for many real networks, the degree distribu- 
tions are different from exactly power laws |8, 9]. For exam- 
ple, the scientific collaboration networks can be better charac- 
terized by the power-law distributions with exponential cutoff 
lllOll . the degree distributions of the email networks llllll . some 
collaboration networks |12(1 . and online user-object bipartite 
networks |13] obey the stretched exponential forms, and the 
double power-law distribution seems a better way to describe 
the air transportation networks 1 14- 16]. In this Letter, we fo- 
cus on the Mandelbrot law or called the shifted power law 
ilTll. which can be written as: 



p{k) (X {k + c)~ 



(3) 



TABLE I: Fitting exponents and errors for power law and Mandelbrot 
law on six real networks. P-error and M-error stand for square errors 
from power-law fitting and Mandelbrot-law fitting, respectively. The 
smaller the error, the better the fitting. 



Networks 7(c = 0) P-error 7 



M-error 



SC-Small 


2.4 


5.69 


4.3 


10.5 


1.81 


SC-Large 


2.3 


6.38 


3.6 


10.3 


2.00 


PPI 


3.0 


1.94 


2.9 


-0.1 


1.90 


Slashdot 


1.7 


3.99 


1.9 


1.2 


3.33 


USAir 


1.6 


3.56 


1.7 


-0.9 


3.29 


UCI 


1.8 


3.74 


1.9 


0.4 


3.63 



where 7 is the power-law exponent and c is the shifting co- 
efficient. In fact, for the well-known BA model, the degree 
distribution 



pik) 



2m{m + 1) 



2m'^k~ 



(4) 
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(fc + 2)(fc-|-l)fc 

obtained by the master equation [18], is not an exactly power- 
law distribution. This distribution can be approximated as 
p{k) oc (fc + 1)^'^, which also satisfies the Mandelbrot law 
with 7 = 3 and c ~ I. 

Recently, the Mandelbrot law has been applied to charac- 
terize the degree distributions of some real networks iflOluOll . 
Here we test the validity of the Mandelbrot law on six real 
networks: (i) SC-Small.-A scientific collaboration network 
according to e-print manuscripts on condense matter physics 
from 1995 to 1999 in arxiv.org [10]; (ii) SC-Large. -Similar 
to SC-Small but based on manuscripts published from 1993 
to 2003 1I21II: (iii) PPI.-A protein-protein interaction network 
of yeast ||22I1 : (iv) Slashdot. -A social network consisted of 
friend/foe links, attracted from the online service website, 
Slashdot 1I23I1 . (v) USAir. -An air transportation network in the 
United States 12411 . (vi) UCI. -A social network of students at 
University of California, Irvine |25]. We try the least square 
method for both power law and Mandelbrot law, and the fitting 
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FIG. 1: Illustration of power-law and Mandelbrot-law fittings for two 
scientific collaboration networks. The left and right plots respec- 
tively display the networks of smaller and larger sizes (i.e., SC-Small 
and SC-Large). Obviously, the Mandelbrot law performs much bet- 
ter. 



exponents as well as square errors are shown in Table 1 . From 
this table, we could conclude that: (i) Using the Mandelbrot 
law can generally improve the fitting accuracy compared with 
the power law since the former has one more coefficient; (ii) 
Sometimes the two fitting methods give more or less the same 
errors, and in these cases, the two power-law exponents are 
almost the same and the shifting coefficient is usually close to 
zero; (iii) Sometimes applying the Mandelbrot law can largely 
improve the fitting accuracy, and then the two power-law ex- 
ponents are far different while the shifting coefficient is much 
larger than zero and its significant role cannot be neglected. 
Figure 1 displays the two cases with remarkable differences 
between two fitting methods, from which the advantage of the 
Mandelbrot law is demonstrated. 

A number of tools have been developed to get the analyt- 
ical solutions of network degree distributions, including the 
mean-field approximation, the master equation, the rate equa- 
tion, and so on |7, 18, 26-28]. Most of these known analyti- 
cal methods only concentrate on the power-law exponent, yet 
pay less attention to the value of shifting coefficient, which, 
however, plays a significant role in determining the shape of 
degree distributions (see figure 1). Even worse, we will show 
later that the widely used difference approximation, no mat- 
ter forward difference or backward difference, will result in 
considerable deviations to the real value of the shifting coeffi- 
cient. 

We here investigate a model embodying a linear preferen- 
tial attachment, which can be considered as an extension of 
the famous BA model. Initially, our model starts with a fully 
connected network with too nodes and mo(mo — l)/2 links. 
If the final network size is S, it should satisfy the condition 
Too <C S. After initialization, at each time step, a new node 
will be added into the network, which will connect to to old 
nodes. The probability of an old node i to be connected is 
linearly correlated with its degree ki, say 



n(fc,)-^(afc^ + /3), 



(5) 



where a and /3 are two parameters, and N is the number of 
nodes at that time step. The self-loop and multiple links are 
not allowed. This model will degenerate to the BA model if 
/? = 0. The parameters a and (3 satisfy the normalization 




FIG. 2: Degree distribution of the modeled network with /? = 
and m — 3. Since /3 = 0, it is equivalent to a BA network. Com- 
pared with the results of backward difference approximation (blue 
dash line) and forward difference approximation (green dot line), 
the present solution (red solid line) is closer to the simulation re- 
sult (round donuts). The network size is 5 = 10 and the results are 
obtained by averaging over 100 independent realizations. 



condition 



^7r(fcMfc)-l, 



(6) 



where 7r(fc) = A^n(fc) and TT{k)p{k) is the probability a se- 
lected old node is of degree k. That is 

J2 7r(fc)p(fc) = ^(afc + /3)p{k) = 2ma + /3 = 1. (7) 



It is equivalent to: 



a 



2^^'-^^- 



(8) 



The rate equation is based on the assumption that the added 
nodes and links, during a time step, have no influence on the 
degree distribution of the network, namely in the thermody- 
namic limit, the degree distribution approaches to a steady 
form. Denoting p{k) the steady degree distribution and N the 
number of nodes in the current time step, if N is large enough, 
then the number of nodes with degree k is approximated to 
Np{k). Analogously, the number of nodes with degree k in 
the next step is {N + l)p{k). With to, links added, the number 
of nodes with degree k in the time step A^ + 1 reads 

{N+l)p{k) = Np{k)+m7r{k~l)p{k-l)~m7r{k)p{k)+Skm, 

(9) 
where imT{k — l)p(fc — 1) and imT{k)p{k) represent, respec- 
tively, the number of nodes whose degree changes from fc — 1 
to k and from fc to fc + 1 in this time step. 6km accounts for 
the specific degree equal to m, namely Skm — 1 when k — m 
and 5km = otherwise. Eq. (9) is the familiar form of the 
well-known rate equation ll26ll27ll . which is usually solved by 



the difference approximation, however, here we will show a 
different analytical method, and will later compare our results 
with the ones obtained by the difference approximation. 
Eq. (9) is equivalent to 



I p(m) — T-; — ^ — T, k = m, 

I p(fc) [1 + r7T,7r(fc)] — rmr{k — l)p{k — 1), fc > m. 

Reminding the linear relation 

TT{k) — ak + /3, 



(10) 



(11) 



considering Eq. (8), the probability of a newly added link 
connecting to an old node with minimum degree is 



Tr{m) 



2m 



m + /3 



1 + /3 



(12) 



Clearly, this probability should be no less than zero and no 
larger than one, and thus —l<f3<l- According to Eq. (10), 
the probabiUty density of m-degree nodes is 



p{m) = 



1 



1 + mn{m) 2 + m(l + /3) ' 
Substituting Eq. (8) into Eq. (10), we get 



(13) 



p{k) 



k 



2(1 + ml3) 



1-/3 



2mf3 
1^ 



1 



Specifying: 



2m/3 



1, b. 



2(1 + to/ 



p{k-l). (14) 



(15) 



1-/3 ' l-;3 ' 

then Eq. (14) can be rewritten in a recursive formula as 

k + a 



p{k) 



k + b 



p{k-l). 



Taking logarithm in both sides of Eq. (16), we get 



log 



p{k) 
p{k-l) 



log 



k -\- a 
k + b' 



(16) 



(17) 



With the ansatz that p{k) follows the Mandelbrot law, sub- 
stituting Eq. (3) into Eq. (17), we obtain the relationship 
between power-law exponent 7 and shifting coefficient c as 



log- 



k + b 
which is equivalent to: 

1 



a k — 1 

= 7 log , , 



log 



1 + H 



af 1 

f = 7 log - 



1 



4 



(18) 



(19) 



Under the approximation with large fc, through the second or- 
der Taylor expansion of Eq. (19) with 1/fc being the variable, 
we can get the power-law exponent 



7 



1 



1-/3' 



(20) 
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FIG. 3: The comparison of degree distributions with different shift- 
ing coefficients given m — 5 and S — 10000. Compared with the 
case of /J = 0.8 and shifting coefficient c = 50.4 (green down- 
triangles), the degree distribution of the none-shifting case with 
P = —0.1 and c = (purple circles) is much closer to a straight 
line in the log-log coordinates. The results are obtained by averaging 
over 100 independent realizations. 



and the shifting coefficient 



b + a + 1 1 + 2to/? 
2 ^ 1-/3 ■ 



(21) 



Eq. (20) and Eq. (21) declare that 7 only depends on /3, 
while c is related to both /3 and m. When to/3 is very large or 
j3 -^ 1, a and b are both very large, and the Taylor expansion 
cannot be applied on Eq. (19). Under such condition, Eq. (16) 
can be approximately rewritten as 



p(fc)«-p(fc_l), 



(22) 



then the degree distribution is close to an exponential form. 
It is easy to be understood since when /3 — )► 1, the selection 
of old nodes is almost random. When B — 0, a ~ tt-, our 
model degenerates to the BA model, and we get a — ~1, 
fe = 2, 7 = 3 and c = 1, with degree distribution being 



pik) 



l/'(2,TO+l) 



(fc + 1) 



-3 



where 



^(x) = T'ix)/T{x) 
is the Digamma function with 



r{x) ^ I e^'e-^dt 



being the Gamma function and 



ip{n,x) 



d"Tp{x) 
dx"- 



(23) 



(24) 



(25) 



(26) 



We next compare the present method with methods based 
on the difference approximation. We first introduce the back- 
ward difference approximation, which assumes 



^=p{k)-p{k-l). 
Substituting Eq. (27) into Eq. (16), we get 



pik) 



which is equivalent to 



k + a 
k + b 



p{k) 



dp 
Ik 



(27) 



(28) 



dp 
dk 



a — b 



p{k) 



k + a 
that leads to the solution 

Similarly, the forward difference approximation assumes 
dp 



(29) 



■£=p{k + l)-p{k), 



and then Eq. (16) can be rewritten as 

k + a + 1 



p(fc + l) = 



1 



p{k), 



-b 



-p{k). 



which is equivalent to 

dp 

dk^ k + l + b' 

In this case, the solution is 

p{k) oc(fc + 6 + l)-(''-"). 



(31) 



(32) 



(33) 



(34) 



The three methods all indicate that the Mandelbrot law will 
emerge from an evolving network with linear preferential at- 
tachment, and give the same power-law exponent j = b — a. 



In contrast, the shifting coefficient are different: c^':^^^'^^ = 

a+l+l^ ^backward ^ ^ ^^^j ^forward ^ fe + 1. !„ Fig. 1, 

we compare the degree distributions obtained by these three 
methods with the simulation results and show that the present 
method is more accurate. 

Although we usually refer to the concept of scale-free net- 
works, neither the BA networks nor most real networks have 
very precisely power-law degree distributions. The present 
method suggests that we can obtain a more precisely power- 
law distribution by setting a right (3 that corresponds to a zero 
shifting coefficient. Since the degree distribution is 



,,, ,, 1 + 2to/3V+t^ 
p(k) oc l^k+ ^_ 



(35) 



(30) the non-shifted degree distribution asks for 



0, 



(36) 
Therefore, given 



l + 2m/3 

namely 13 — ~^ and p{k) oc k 
the linear preferential attachment, the non-shifted power-law 
exponent is determined by the average degree and can never 
exceed 3. Figure 2 compares two degree distributions, respec- 
tively with c — and c = 50.4, from which one can con- 
firm that the non-shifted power law is indeed much closer to a 
straight line in the log-log coordinates, and the shifting coef- 
ficient largely affects the shape of degree distribution. 

In summary, we extend the BA model to an evolving model 
with linear preferential attachment and show that the corre- 
sponding degree distribution obeys the Mandelbrot law. The 
shifting coefficient, usually being ignored in the literature, 
largely affects the shape of degree distribution. In puzzlement, 
the backward and forward difference approximations will lead 
to different solutions on shifting coefficient. Our analysis in- 
dicate that both of them are inaccurate, and we propose an 
analytical method that results in a more accurate solution. 
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